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Enhancement and suppression of nonlinear processes in coupled systems of plasmonic converters 
and quantum emitters are well-studied theoretically, numerically and experimentally, in the past 
decade. Here, in difference, we explicitly demonstrate -with a single equation- how the presence of 
a Fano resonance leads to cancellation of nonresonant terms in a four-wave mixing process. Cancel¬ 
lation in the denominator gives rise to enhancement in the nonlinearity. The explicit demonstration, 
we present here, guides us to the method for achieving more and more orders of magnitude enhance¬ 
ment factors via path interference effects. We also study the coupled system of a plasmonic converter 
with two quantum emitters. We show that the potential for the enhancement increases dramatically 
due to better cancellation of the terms in the denominator. 
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Introduction. Plasmonic resonators, such as metal 
nanoparticles (MNPs), graphene nano-islands [T] and 
transparent conducting oxides [2], confine the incident 
optical electromagnetic field into nano-dimensions in the 
form of surface plasmon (SP) excitations. When a quan¬ 
tum emitter (QE) is placed on one of these hot-spots, 
plasmons interact with the emitter strongly [3]. Due to 
the small decay rate of the QE, Eano resonances appear 
in the spectrum of the plasmonic material, where absorp¬ 
tion vanishes. Eano resonances, appearing in the linear 
response [4], can extend the lifetime of plasmons [5] and 
lead to further enhancement of the local field [6] which 
makes coherent plasmon emission (spaser) possible [7]. 

Linear Eano resonances do not emerge only when a 
plasmon is coupled with a QE, but also observed for two 
interacting plasmonic materials if one has a smaller decay 
rate [8]. The origin of Eano resonance is analogical to 
electromagnetically induced transparency (EIT) [9], and 
can clearly be understood in terms of interference of two 
absorption paths m- 

In the past decade, the role of Eano resonances in the 
modification of nonlinear response is also well-studied. 
Presence of a quantum emitter, with an energy spac¬ 
ing uJeg close to the output frequency, can enhance or 
silence the frequency conversion [nHU]. Both effects, 
suppression and enhancement, are desirable for devices 
designed to operate in the linear and nonlinear regimes, 
respectively. Operation of high-power lasers and fiber op¬ 
tic cables [18] necessitates the avoidance of second har¬ 
monic generation (SHG) and Raman scattering, in order 
to prevent the loss of energy to other modes. On the 
other hand, nonlinear imaging with SHG m, four-wave 
mixing (EWM) [20| and Raman spectroscopy [21], and 
generation of nonclassical/entangled plasmons [22] desire 
enhanced nonlinear response. 


Eour-wave mixing (EWM) process -beyond its tech¬ 
nical applications such as superresolution imaging m, 
ultrafast optical switching [23] and nonlinear negative 
refraction [24]- is also important for studying one of 
the fundamental inquiries in quantum information. In 
a EWM process, two fields become entangled in addition 
to squeezing in one of the fields. Hence, EWM is a spe¬ 
cial process in which the interplay between single-mode 
nonclassicality and two-mode entanglement coexist nat¬ 
urally [25]. Localization is shown to enhance also EWM 
in plasmonic materials [26] due to the larger overlap in 
mode-integrals for that process [laiiH]. 

Similar to SHG [13] , presence of quantum emitters at¬ 
tached to plasmonic converters modifies the EWM pro¬ 
cess [29H33] . Analogous to linear response [8] and SHG 
[in [HUH HU, coupling of two plasmonic materials with 
different decay rates creates modifications due to Eano 
resonances also in the EWM process [34]. Eano reso¬ 
nances in EWM can show themselves even as the coupling 
between superradiant and subradiant collective modes of 
plasmonic nanoclusters [33], similar to dark states in a 
plasmonic SHG process [35] . 

Presence and behavior of such resonances are well- 
studied in the literature. In this Letter, in difference, 
we demonstrate the path interference effects explicitly in 
a plasmonic EWM process. We provide a single equation 
[see Eq. ( [Tq] )] for the steady state amplitude of the EWM 
conversion. Interaction with the quantum emitter intro¬ 
duces extra terms in the denominator of the conversion 
amplitude. We show that enhancement emerges simply 
due to the cancellation of these extra terms with the non¬ 
resonant term in the denominator. On the contrary, con¬ 
version suppression emerges when the extra term grows 
several orders of magnitude -due to the small quantum 
decay rate- and makes the denominator blow up. 
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FIG. 1: (Color online) (a) A quantum emitter (QE) is placed 
to the hot-spot of a plasmonic converter. Field localization 
provides a strong plasmon-QE interaction, (b) Two lasers of 
frequencies uji and UJ 2 drive the plasmon modes of resonances 
and Q 2 , respectively. 2 plasmons in the mode, oscil¬ 
lating with cji, combine and decay into two plasmons with 
different frequencies, one into Q 2 mode and another one into 
the Q 3 mode. Decay into the Q 2 mode is stimulatd by the UJ 2 
laser. Frequencies are scaled with uji. 


We also study the system where a plasmonic converter 
is coupled to two quantum emitters (QEs). We show 
that a better cancellation of the nonresonant term -this 
time also the decay term- can be achieved in the FWM 
process, see Eq. ( p!4| ). FWM enhancement increases by an 
order of magnitude (15 times) compared to the system 
where the plasmonic converter is coupled with a single 
QE. 


We present the second quantized Hamiltonian for the 
system of a plasmonic converter (FWM) coupled to a sin¬ 
gle quantum emitter. We obtain the equations of motion, 
introduce decay terms, convert operators to c-numbers 
(since not interested in entanglement) and obtain the 
steady-state result for the converted FWM frequency [see 
Eq. (10)]. We repeat the similar procedure for the sys¬ 
tem of plasmonic converter coupled to two QEs. In a 
previous publication we showed that one does not need 
to concern about the retardation effects. We simulated 
SHG process in a coupled system using 3D Maxwell equa¬ 
tions and we observed that the spectral position and the 
strength of Fano enhancement are not affected signifi¬ 
cantly (see Fig. 5 in Ref. ini)- 


Hamiltonian. Nonlinear processes in plasmonic con¬ 
verters take place through plasmons [36] due to growing 
mode (overlap) integrals [27l [28] with the localization. 
In FWM process, + H.c.^ the mode integral de- 

terming the strength of this process is proportional as 


m 


Xpwm J dh E;{r)E;{r)E^,{r) , (1) 

where Ei are the spatial extents of the electric (or po¬ 
larization) fields of the plasmon modes around the 
nanoparticle. A nonvanishing Xfwm necessitates some 
parity requirements for the fields of the plasmon modes, 
in symmetric particles. 

The dynamics of the conversion process can be de¬ 
scribes as follows. Plasmonic converter has three reso¬ 
nances, of frequencies , and in the relevant 

region of the spectrum (see Fig. IT)- Two lasers of fre¬ 
quencies uji and UJ 2 drive the two plasmon modes 
respectively. The lasers excite surface plasmons in the 
di-mode (Qi) and d 2 -mode {^ 2 )^ wherein plasmons os¬ 
cillate at uJi and co’ 2 . 2 plasmons in the fti mode combine 
and decay into two plasmons with different frequencies, 
one in the Q 2 mode and another one in the Qs mode. 
The plasmon generated in the ds mode (f^s) oscillates 
with frequency 003 = 2uJi — Co’ 2 - Decay into Q 2 mode is 
stimulated by the presence of plasmons in this mode. Co’s 
plasmon oscillations emerge due to the energy conserva¬ 
tion and the presence of a plasmon mode at fts ^ Co’s. In 
order to observe the Co’s photons in the far field, plasmons 
in the Q 3 mode can be transformed to photons selectively 
by gratings m- 

Therefore, the Hamiltonian for the total system can be 
written as the sum of the energies of the quantum emitter 
(Hq)^ plasmon modes (^pis), the plasmon-QE interaction 

(^int) 


Ho = hu)e \e) (e| + hwg \g) {g\ , (2) 

Hp\s = hQ.id\di + hQ.2d\d2 + hQ,sa\ds , ( 3 ) 

Hint = fi(^fdl\g) {e\+f*d3 \e) {g\^ , (4) 

and including the two laser pumps (^^pump) and the non¬ 
linear FWM process (^^fwm) 

Hp = ih (a\ + 4 + iJ.c.) , (5) 

-^FWM ~ ^Xfwm (^3^2^i 4“ ) 5 (6) 

where / is the strength for plasmon-QE interaction, uje 
(ujg) is the excited (ground) state of the QE, 5p^^ and 
are the strengths of the two lasers of frequency uji and 
Co’ 2 , and mode-integral of the nonlinear process Xfwm is 
defined in Eq. 0 . We choose the energy level spacing of 
the QE, ujeg = uJe — ojg^ such that it falls into the spectral 
region of the plasmon mode. Hence, we consider the 
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interaction of QE with as mode only, which provides sub¬ 
stantial simplification making analytic results [Eqs. ( [1Q| ) 
and (14)] possible [38] . 

Equations of Motion. Dynamics of the plasmon opera¬ 
tors CLi and the density operators for the QE peg = \e){g\ 
and pee = k)(^l can be obtained from the Heisen¬ 
berg equation of motion, e.g. ihoi = [a^,^]. We are 
not interested in the correlations in the system. Oscilla¬ 
tion amplitudes (or occupation numbers) of the plasmon 
modes can as well be handled by replacing operators by 
c-numbers, i.e. Oi ai. When we introduce the plas- 
monic linewidths 71 , 72 , 73 and the decay rate of the 
quantum emitter jegi dynamics is governed by coupled 
equations 

di = - 71 ) Qfi - i2xFWMC^lc^2C^3 + (7a) 


d 2 = (-ifl 2 - 72 ) 0^2 - (yp) 

d 3 = (-ifl3 - 73 ) <^3 - iXFWM(^2^l - ifpge (7c) 


Pge — ( i^eg leg) Pge '^/<^3 {Pee Pgg) (7^) 


Pee — XeePee {f^sPge f ^SPeg) (7c) 


with pgg = 1- Pee and Xee = 27 , 


eg- 


Steady-state. When the dynamics of Eqs. (7a)-(7e) is 


examined by placing exponential solutions, one identifies 


(0 — 


—iujit 


and 


Pge Pge^ 


—i(2uji —uj2)t 


(8) 


as the steady-state oscillations, where and pge are con¬ 
stants. It is worth noting that, in order to be confident 
about the validity of Eq. j® , we obtain Eigs. 2 and 3 by 
time evolving Eqs. (7a)-(|7e|). When we place Eq. ^ into 


Eqs. (7a)-(7e), we reach the coupled nonlinear equations 


[i {fti - uji) + 7 i] di + i2xFWMd^d2d3 = 4^^ (9a) 

[i {^2 - ^ 2 ) + 72 ] d 2 + iXFWMdgdi = 4^^ (9b) 

[i {fts Euj 2 - 2a;i) + 73 ] d3+ixFWMd2di = -ifpge (9c) 
[i {p^eg H“ ^2 2cJi) -|- Xeg\ Pge — '^f^S {Pee Pgg) (9d) 

leePee = i {f^lPge “ F^Speg) (9e) 

for the steady-state plasmon occupations, jd^p 

The single equation. With a simple algebraic manipu¬ 
lation, using Eqs. ( [^ and ( |9d[ ), one can reach the very 
simple and useful equation 

iXFWMd2d^ 

^3 I P|2 

[iK,+J727)+7e.] - (^3 + ^^2 - 2 w ,) + 73] 

( 10 ) 




FIG. 2: (a) Relative enhancement of the FWM process in 

the presence of coupling to the QE, compared to bare plas- 
monic converter. The firs term in the denominator of Eq. (10) 
cancels the following nonresonant term. Such a cancellation 
occurs for uoeg — 1-53 [Eq. which yields an 80 times en¬ 
hancement. (b) When ujeg = 2.uji —uj 2 = the extra term in 
the denominator becomes y\ff‘/xeg very large and suppresses 
FWM. 


which determines the number of FWM plasmons |d 3 p, 
with y = pee — Pgg IS the population inversion for the QE. 
When coupling between the plasmonic converter and the 
QE is not present, / = 0, Eq. (10) displays the simple 
resonance condition for EWM. Maximum conversion is 
attained when 20 ;! — UJ 2 = ^ 3 , as should be expected. 
This happens when the converted frequency is resonant 
with the plasmon mode. In Eq. ( [Tq| ) we assume resonant 
pumping for the two lasers, uji = Oi and UJ 2 = for 
simplicity. We refer the term — (20;! — 002)1 present in 
the denominator of (10), as the nonresonant term. EWM 
is a very weak process. Hence, one can safely consider 
that occupations of the two plasmon modes, |dip and 
|d 2 p, are almost unaffected from the conversion process. 
In most cases, population stays weakly excited and one 
can imagine it as ^ — 1. 

Enhaneement in EWM, in |d 3 p, can be achieved if 
one arranges the imaginary part of the first term of the 
denominator, \f^ y/[i {ooeg P- 0 J 2 — 2oji) -h Xeg\i fo cancel 
the nonresonant term (fl 3 002 — 2uji). Therefore, if the 
level spacing of the QE is chosen such that 


OJ 


* 


W 3 + 


\f?y 

^3 ~ ^3 


l/IV 


(^3 — ^ 3 )^ 


— 47 ^ 

/e^ 


( 11 ) 


there emerges an enhancement peak in the EWM spec¬ 
trum, where 0)3 = 2uji — 002 . This is depicted in Eig. [^. 
Erequencies are scaled with ooi and cji = Oi = 1, 
0 J 2 — f ^2 ~ 0.5, oj^ = 1.5, O 3 = 1.85, Xfwm ~ 10 
/ = 0.1, 7l = 72 = 73 = 0.01, Jeg = 10“®. 
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Suppression of FWM, on the contrary, can be obtained 
by arranging the term \ ff y/[i {ujeg + - 2 cJi) + ^eg] to 

blow up. This happens when the level spacing of the QE 
is set to uJeg = 2co’i — co’ 2 . Since ^eg ^ lO^Hz and plas- 
mon resonances are typically ^ lO^^Hz, the extra term 
becomes very very large and leads to vanishing FWM, 
jdsp. This is depicted in Fig. [^. 

More and more enhaneement. The enhancement in 
Eq. ( [1Q| ) is limited by the presence of 73 , which is not a 
small quantity. Hence, one naturally seeks a way to in¬ 
troduce more extra terms to get rid of (cancel) the decay 
term 73 in the denominator. For this purpose, we con¬ 
sider the system where the plasmonic converter is coupled 


to two QEs, with strengths fi and / 2 . The two QEs also 
interact with each other with a strength g. 

Hamiltonian needs to be modified with addi¬ 

tional and replacing terms 


ffqE-QE = fi-( 5 |e 2 ) ( 52 I ® I 51 ) (eil + H.C.) , (12) 


ifint = 


fi (/i4 l^i) (eil + 724 152 ) ( 62 ! + ifi.c.) . (13) 


When the procedure, similar to Eqs. (7a)-(7e) and (9a)- 


, is applied m one again reach a single equation 


0^3 = 


^Afwm yh'ih'2 -r yiy2 \y\ j 

{yi \fif + y2 1/2!^ ^1) + iym ififig* + A*As) - 6 (/^lA + ym Isl^) 


a2«i ) 


(14) 


where ^3 = i {fts ^ UJ 2 —+ 73 and /3j = 

i + CJ 2 — 2co’i^ with j = 1,2 enumerates QEs. 

(j) rj\ 

Here, uJeg and ^eg are the energy level spacing and the 
decay rate for the QE. yi and ^2 are population in¬ 
versions. 


This time, the extra terms cancel 73 more efficiently. 
In Fig. we see that the cancellation in the denomina¬ 
tor [Eq. results 1200 time enhancement for FWM, 
which is 15 times larger compared to the one for single 
QE [Eq. [^. In order to maximize 0^3 in Eq. (14), we 
numerically minimize the denominator by varying /i, / 2 , 
^5 ^eg and id eg m- In principle, the desired phase for 
the QE-QE interaction (g) or plasmon-QE interactions 
if 1 , 2 ) can be obtained by arranging the phase of the in¬ 
teraction integral [28] by reshaping the two QEs. In the 
case of all-plasmonic Fano resonances naiii] such an ar¬ 
rangement among the interacting plasmon modes would 
be much simpler. 


Conelusion. We obtain the steady-state amplitude of 
FWM process as a single equation for a plasmonic con¬ 
verter coupled to quantum emitters (QEs). The denomi¬ 
nator of this equation reveals the path interference effects 
explicitly. Presence of QEs, interacting with the plas¬ 
monic converter, introduces extra terms in the denomi¬ 
nator. When the extra terms cancel nonresonant term 
for the conversion, enhancement of FWM is observed. 
Contrary, when the extra term blows up due to the long 
lifetime of the QE, suppression of FWM is observed. 


We utilize this observation as a tool for obtaining larger 
and larger enhancement in FWM. We show that, for the 
coupling of the plasmonic converter to two QEs, a better 
cancellation of the denominator is achieved. Enhance- 
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FIG. 3: (Color online) When the plasmonic converter is cou¬ 
pled to 2 QEs, there introduces additional terms, see Eq. (14). 
These terms manage to partially cancel also the 73 term, 
which is not possible using a single QE as in Eq. (10). En¬ 
hancement reaches to 1200 times. 


ment increases by a factor of 15 on top of the resonance 
obtained by coupling with a single QE. Fano resonances, 
both in the linear and nonlinear response, emerge also by 
coupling the plasmonic converter with a plasmonic ma¬ 
terial which has a smaller decay rate [laiiiHisi. Hence, 
our method is possible to be generalized [39] to interact¬ 
ing plasmonic clusters, since the experiments are easier 
to conduct with nanoparticles. 

M.E.T and S.K.S. acknowledge support from 
TUBiTAK-1001 Grant No. 114F170. 
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